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Abstract 

We study the magnetoconductivity induced by the axial anomaly via the 
chiral magnetic effect in strongly coupled holographic models. An important 
ingredient in our models is that the axial charge is non-conserved beyond the 
axial anomaly. We achieve this either by explicit symmetry breaking via a 
non-vanishing non-normalisable mode of an axially charged scalar or using a 
Stiickelberg field to make the AdS-bulk gauge field massive. The DC magneto- 
conductivites can be calculated analytically. They take a universal form in terms 
of gauge field mass at the horizon and quadratic dependence on the magnetic 
field. The axial charge relaxation time grows linearly with magnetic field in the 
large B regime. Most strikingly positive magnetoconductivity is still present 
even when the relaxation times are short rs ~ l/(7rT) and the axial charge can 
not be thought of as an approximate symmetry. In the U{1)a explicit breaking 
model, we also observe that the chiral separation conductivity and the axial mag¬ 
netic conductivity for the consistent axial current vanish in the limit of strong 
symmetry breaking. 
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1 Introduction 


Anomaly induced transport phenomena have been in the focus of much theoretical and 
experimental research lately. One particular example is the so-called chiral magnetic 
effect (CME) [1].^ It states that in the presence of an imbalance in the number of 
left-handed and right-handed charged fermions an electric current is generated in a 
magnetic held whose origin can be traced back to the presence of the axial anomaly 

d.JS = ( 1 . 1 ) 

with the anomaly constant c = l/(27r^) for a single Weyl fermion. The chiral magnetic 
effect is most conveniently expressed in the form 

J = c/i5B, (1.2) 

where the axial chemical potential ps parametrises the imbalance in the number of left- 
and right-handed fermions. The subtleties in the interpretation of this formula have 
been discussed extensively in the literature [5, 6, 7]. Here it shall suffice to emphasise 
that ps is a parameter that characterises a state in which the chiral imbalance has 
been induced in a dynamical way, e.g. via the very axial anomaly by switching on 
parallel electric and magnetic helds. Indeed such a field configuration will induce an 
axial charge whose Fourier transform is 

Ps = —cE ■ B . (1.3) 

UJ 

We also relate axial charge and axial chemical potential via an axial charge suscepti¬ 
bility /is = xsps and using (1.2) we obtain the electric held response in the current 

7 r'^B B 

J^ = aEE, + -^^E^. (1.4) 

^ X5 

Here cte denotes the quantum critical conductivity (also named as charge-conjugation 
symmetric conductivity, or Ohmic conductivity). In [8] this effect has been studied 
in both hydrodynamics and holography and the result indeed conhrmed the above 
formula. It was found that the formula applies even in the large magnetic held regime in 
which the anomalous hydrodynamics of [3] is not applicable anymore. More precisely it 
turned out that cte became anisotropic and its longitudinal component is exponentially 
suppressed for large B, see (A.2). At the same time the axial susceptibility becomes 
linear in the magnetic held, consistent with the expectation of Landau level physics 

^Early version of the CME have appeared in [2]. It has been derived in a variety of approaches 
ranging from hydrodynamics to the gauge gravity duality, e.g. [3, 4]. 
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at weak coupling. The delta function peak indicated by the imaginary pole in (1.4) 
therefore showed a cross over to a scaling behaviour that was linear in B. This indeed 
is consistent with weak coupling results [9, 10, 11]. 

Thus if axial charge was conserved up to the axial anomaly this would result in a 
perfect superconductor. The pole in the imaginary part implies indeed a delta function 
singularity located at zero frequency in the real part via the Kramers-Kronig relation. 

Alas nature seems to abhor axial symmetries even beyond the electromagnetic axial 
anomaly. Fermions are typically massive. Another way the axial symmetry is broken 
is via the QCD contribution to the anomaly.® In both cases it is more realistic to allow 
for a non-conservation of the axial charge and introduce an axial charge relaxation time 
Ts. In this case the frequency is shifted a; —>■ a; + i/rs in (1.4) resulting in a hnite DC 
magnetoconductivity 

Ji = aEEi + T5 ^ Ej . (1.5) 

X5 

In the studies up to now the relaxation time has been introduced via a naive 
relaxation time approximation. Basically it was added “by hand” and did not represent 
dynamics inherent to the holographic system. 

Recently a number of experiments on Dirac or Weyl-(semi-)metals have reported 
strong positive magnetoconductivity along the magnetic held [13, 14, 15, 16, 17]. They 
furthermore always show a quadratic dependence on the magnetic held strength.^ 
These materials are candidates for realising the axial anomaly in a solid state setup. 
They do however also have an intrinsic axial charge relaxation time since the axial 
charge is realised as an accidental low energy symmetry. 

This motivates us to study the anomaly induced magnetoconductivity in models in 
which the relaxation time is an intrinsic property of the model. In particular we will 
concentrate on two models with diherent symmetry breaking mechanisms. 

The hrst model uses a massless bulk gauge held for the axial current. We break 
the axial symmetry by introducing a tree level coupling to a non-normalisable mode 
of an axially charged scalar held. We chose the bulk mass of the scalar so that the 
dimension of the breaking parameter is one and can be interpreted as a mass M for 
the dual fermions. The weak coupling interpretation is now that the divergence of the 

®We assume the electromagnetic fields to be external not dynamical contributions in the axial 
anomaly. Instabilities arising from dynamical gauge field have been discussed in [12]. 

^ There is a regime for very small magnetic field in which the conductivity first decreases with 
magnetic field. This seems to be a disorder effect called “weak anti-localization” and most likely 
represents physics unrelated to the anomaly. 
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axial current is of the form 

d,J^ = M^75V' + . (1.6) 

Although the system has the same ingredients as a holographic superconductor its 
dynamics is quite different due to the fact that it is the non-normalisable mode that 
is non-vanishing. In particular the DC conductivity is hnite and the quasinormal 
mode spectrum does not contain sound modes. The axial charge relaxation time is 
determined by the the mass parameter M. When t^T 3> 1 we can think of the axial 
symmetry as approximately conserved. In this regime we indeed hnd linear scaling of 
the susceptibility and at the same time linear scaling of the relaxation time Ts with large 
magnetic held. Consequently the magnetoconductivity scales quadratic with magnetic 
held which is at odds with the naive relaxation time approximation that assumes 
independence of on B. When t^T ~ 1 we do not really have justihcation anymore 
to think of the axial charge as approximately conserved. Generically the quasinormal 
frequencies for non-conserved operators have imaginary values of the order of T and 
therefore there is no hierarchy between the relaxation of the axial charge and a generic 
perturbation. In this regime the derivation of the magnetoconductivity (1.5) outlined 
above becomes invalid. Somewhat surprisingly we still hnd a magnetoconductivity 
that depends quadratically in the magnetic held strength in this regime! In fact we 
can hnd an exact analytic expression for the magnetoconductivity that conhrms the 
exact quadratic dependence on the magnetic held for all values of the charge relaxation 
time. 

The second model uses a massive gauge held. As is well-known global symmetries 
correspond to massless gauge held in holography. The mass of the gauge held gives the 
dual current operator an anomalous dimension A = — 1 -|- \/l + m?. Since conformal 
symmetry implies that the dimension of a conserved current is precisely d in d + 1 
dimensions this means that a massive gauge held corresponds to a non-conserved vector 
operator in a conformal theory. In QCD the axial anomaly for U{1)a receives not only 
electromagnetic but also a gluonic contribution. Holography ehectively replaces the 
gluon dynamics by gravity in anti-de Sitter space. It has been argued some time ago in 
[18] that the gravity dual of this gluonic contribution to the anomaly is a Stiickelberg 
mass term for a vector held. Therefore the divergence of the axial current is dual to 
the weak coupling form 

{c,^F^,Fpp + Cstrongtr ) (1.7) 

where the anomaly coefficients Cem and Cstrong depend on the fermion spectrum.® An- 

®In particular to make sure that the gluonic contribution is not subleading in the supergravity limit 
one needs to assume typically Veneziano limit with the number of flavor to be of the same order as 
the number of colours [19]. 
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other way of viewing such models is by noting that the anomaly is a dimension four 
operator and therefore should couple to a marginal scalar held, this is the Stiickelberg 
held. One might also include another real scalar held that can serve as the dilaton, 
the source for the kinetic term of the non-abelian gauge helds, which also would allow 
to break the underlying conformal symmetry via a dilaton how. Such models have 
already been introduced and studied in [19, 20]. We will rehne that analysis of the 
magnetoconductivity of [20]. We hnd an analytic result for the magnetoconductivity 
and it turns out that it is of the same form as for the previous model, except that 
now the mass of the gauge held is constant throughout the bulk spacetime. Again the 
magnetic held dependence is quadratic independently of the value of the gauge held 
mass. We also calculate the magnetic held dependence of the axial relaxation time and 
the axial charge susceptibility. Whereas the relaxation time shows linear dependence 
in B the susceptibility has a scaling that is determined by the dimension of the ax¬ 
ial current This is an indication that for this model the simple hydrodynamic 

reasoning does not apply straightforwardly since it would predict a dependence of the 
form Ts/xs. We also emphasise that when the current operator has dimension 3 -|- A 
the source for the time component is not an axial chemical potential but rather a 
true coupling in the theory. 

This paper is organised as follows: in section 2 we study the model with an axially 
charged scalar field. We switch on a non-normalisable mode that serves as an axial 
symmetry breaking parameter (fermion mass). We compute the magnetoconductivity 
analytically and check the expression numerically. We numerically study the large B 
behaviour of the relaxation time and the axial susceptibility and hnd linear scaling for 
both in the regime Tts 3> 1. This is the regime in which the hydrodynamic reasoning 
applies and it indicates indeed quadratic dependence of the magnetoconductivity on 
B. In the regime Trs ~ 1 hydrodynamics does not apply anymore since the axial 
charge is not even approximately conserved but the analytic formula is of course still 
valid in this regime. In section 3 we study the model with a Stiickelberg axion. We 
hnd the analytic result for the magnetoconductivity and check it numerically. In this 
model we can also hnd an analytic result for the axial charge susceptibility. It shows 
scaling with B^^^ for large B. Numerically we also hnd linear scaling in B of the axial 
charge relaxation time. In section 4 we present a comparison of the two models and an 
alternative derivation of the DC conductivity based on near far matching method. We 
hnish in section 5 with a discussion of our results. Technical details of the calculations 
are summarised in the appendix. 


4 


2 Holographic U{1)a explicit breaking model 


It has long been known that for a chiral anomalous fluid in presence of magnetic 
held, a large longitudinal magnetoconductivity due to the chiral anomaly is induced. 
[9, 10, 11]. In a previous study in the framework of hydrodynamics [8], in order to 
get a hnite longitudinal DC magnetoconductivity, it was necessary to include energy, 
momentum and also (axial) charge dissipation. However, in the limit of zero densities 
(i.e. fx/{e + p) = + p) = 0), axial charge dissipation alone suffices to give a hnite 

value. 

In this section, we consider the holographic axial charge dissipation ehect in the 
magnetoconductivity by breaking the axial charge U{1)a symmetry with a scalar 
source. We will also introduce two U{1) symmetries and because the electric U{l)v 
symmetry is still conserved, the scalar held only couples to the axial gauge held in the 
bulk. In this way the anomaly caused buildup of axial charge will be compensated 
by the dissipation sourced by this scalar operator resulting in a hnite DC longitudinal 
magnetoconductivity. We will consider the following action® 


S' = y 




1 1 rv 


GT 




( 2 . 1 ) 


with 

- iqA^ 

where the gauge helds V)i and A^ correspond to the vector and axial U{1) currents 
respectively and d) is a complex scalar held with mass rus- Similar models have been 
studied before in e.g. [21, 22] to describe the dual anomalous superconductor with 
U{l)v X U{1)a symmetry where a charged scalar held is introduced to spontaneously 
break U{l)v while not 17(1)^. In this paper we shall turn on a non-zero source as¬ 
sociated to the dual scalar operator in order to break the U{1)a symmetry explicitly, 
which introduces an axial charge dissipation mechanism in our system. 

Note that we do have a conserved current associated to U{l)v-i which means that 
the electric charge is always conserved. With the mass term for the scalar held, the 
dual scalar operator has a scaling dimension A$ = 2 ± and to make sure that 

the scaling dimension of the axial current does not change, has to be negative and 
above the BF bound (i.e. —4 < < 0). Without loss of generality, in the following 

we will choose = —3, which is the most appropriate value as at this mass the scalar 

®We have set the curvature scale L = 1. 
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operator reproduces exactly the scaling dimension of the mass term for free fermionsd'^ 
Thus the conformal dimension of the dual scalar operator is 3 and the corresponding 
source is of dimension 1. This reminds us to the four dimensional free massive fermion 
systems. 

We will study the system in the probe limit, which means only axial charge dissi¬ 
pation is required to get a hnite longitudinal DC magnetoconductivity. In the probe 
limit, fields live in the Schwartzschild black hole background in the bulk 


= r^( - f{r)dt^ + dx^) + 

r^j[r 


f(r) 



( 2 . 2 ) 


with the dual thermodynamical quantities 

e = 3rQ , s = dvrrQ , T = — . 

71 

The equations of motion for these matter fields in the Schwarzschild-AdS's background 
are 


= 0 , (2.3) 

{FrpFp^ + Fr^Fp,) + = 0 , (2.4) 

DpDP^ - Q . (2.5) 

Let us briefly comment on the Ward identity for the current. The dual consistent 
currents are obtained as the variation of the total action with respect to the gauge 
fields. 


JP = hm + Aae^Pf^P^ApFp^) + c.t., (2.6) 

J,= ym + c.t.. (2.7) 

From the on-shell condition, we have 

= 0 . ( 2 . 8 ) 
^ + 3F,fF^) - iq [>l>(B-<*)* - >1>*(B'<*)]) + c.t. 

(2.9) 

course the properties of the holographic system could be far away from the free massive fermion 
system. However it is very intriguing that we do find some similar properties of our holographic model 
compared to free Fermi gas picture. This is similar to the fermionic holographic case with special mass 
of the probe fermion such that the dimension of the dual fermionic operator approaches that of a free 
fermion where the dual system has Fermi surface but not exactly Landau’s Fermi liquid theory [23, 24]. 
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The counterterm contribution is not explicitly shown because it does not add any 
valuable information. The last term in (2.9) contributes to the 1-point function only 
if the non-normalisable mode of the scalar held is switched on. It is tempting to 
interpret it as the contribution of the fermion mass to the Ward identity via 
The covariant version the current is obtained by removing the Chern-Simons term from 
(2.6 - 2.7), i.e. 

= lim + c.t., (2.10) 

= lim + c.t.. (2.11) 

Since the covariant current is the one which has been widely used in the framework 
of hydrodynamics, we will use the covariant current in most parts of our paper and 
briehy comment on the behaviour of the consistent current in subsection 2.3. 

2.1 Background in the probe limit 

Since we are going to study the magneto response, we turn on a background mag¬ 
netic held in the U{l)v sector. We consider the following background with non zero 
components: 

V';. = (V)(r),52/,0,14(r),0) , = (A(r), 0, 0, ^^(r), O) , <|>(r) = 0(r). (2.12) 

It is useful to note that this ansatz is invariant under transformation (r, tq, Vj, W, ^z) 
b~^{r,ro,Vt,Vz, At, A;^), {t,x,y,z) —)• h(t,x,y,z) and B —)■ b~‘^B. We can set tq to be 
1 by choosing b = tq. Plugging this ansatz into the equations (2.3 - 2.5) we obtain 
hve second order ODEs for hve unknown real functions Vt,Vz, At, Az,(f). The equa¬ 
tions can be found in (A.4 - A.8) of the appendix A.2. We can furthermore reduce 
them into four second order ODEs as (A.9), (A.5), (A.6) and (A.8) for four functions 
Vt,(j), At, Az- Once we have the solutions for these four helds, W is totally hxed by 
imposing normalisable boundary condition. 

For = — 3, near the conformal boundary r —)■ cx) we have 




(2.13) 

, V /, Ai ,■ 

— 1-1 —-mr 

T \ 

) +-0 

(2.14) 

Ai ~ /i 5 (^1 -h ^ In r j 

P5 

2r2 ’ 

(2.15) 

A^ ~ So ^ In r j 


(2.16) 


with Ai = |(/r 5 — Sq)^^, pi = Si = —The interpretation for the boundary value is 
the following: y, are the chemical potential and axial chemical potential respectively. 
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M is the source of the scalar operator which explicitly breaks the U{1)a symmetry and 
the corresponding response can be obtained to he O = The parameter 

So = 0 is the analogue of a superfluid velocity. In our context it can also be interpreted 
as the holographic analogue of the separation of the Weyl nodes in momentum space 
of a Weyl semi-metal. 

For this system, the solutions of four second order ODEs (A.9), (A.5), (A. 6 ) and 
(A. 8 ) are specified by eight integration constants, which correspond to the eight in¬ 
dependent parameters at the conformal boundary p, p^, M,ip, sq, S 2 . Near the 
horizon, there are four independent parameters Aj(ro), I^'(ro), A^(ro), 0 (ro) after spec¬ 
ifying regular horizon boundary conditions for At, Vt, 0 and A^. Thus for fixed values 
of m^, q, we have a /ot/r-parameter family of black hole solutions, which correspond to 

P, /i5, M, Sg. 

For our purpose of breaking the axial charge conservation symmetry explicitly at the 
boundary, we focus on solutions with nonzero M, which sources this explicit breaking of 
axial charge conservation. Solutions with hnite chemical potential p and axial chemical 
potential p^ can only be obtained numerically. In fact, as we already know from [ 8 ], 
the inhnite axial charge transfer due to chiral anomaly exists even at both p = 0 
and /is = 0. Thus in this paper we will mainly focus on the p = p^ = 0 case and 
sometimes also the /i = 0, /is 7 ^ 0 and /is = 0, /i 7 ^ 0 cases for simplicity. We also 
set the parameter sg = 0. In the following we list the bulk background solutions for 
/i = /is = 0 , /i = 0 , /is 7^ 0 and /is = 0 , /i 7^ 0 respectively. 


• The solution for /i = /is = 0 can be obtained easily by choosing Vt = Vz = At = 
Az = 0 and in this case we can still have non trivial solution for 0 with source M. 
0 can be solved to be 


0 = 


M 

Y 



'^ Vor[3/4] 

rF[l/4] 


EllipticK 



(2.18) 


where the function EllipticK gives the complete elliptic integral of the hrst kind. 
As /i = /is = 0, the response of the scalar operator is completely determined by 
the source M and we have 


M 

00 


r[i/4] 

\/^F[3/4] 


ro 


~ 1.18rg, 


O 

M 


II3/4F . 

4F[5/4]2 ° 


-0.456r2, 


^^Note that S'ren = S + Sc.t where the counterterm for this special mass is 


- |<i>|2 + l{\ogr^- 1) [\f^ + + \D^<t>\^]) 


and is the induced metric near the boundary. 


(2.19) 


(2.17) 
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with M, O the source and response associated to the scalar operator and 0o the 
horizon value of the scalar held. In this case note that for M = 0 we only have 
the trivial solution 0 = 0. 

The second simple while interesting solution is for Vt = = 0 while At,Vz ^ 0 

(i.e. /i = 0, /xs 7 ^ 0). In this case solutions are governed by the equations (A.9) and 
(A.8) and we can solve the system numerically. The plot for the scalar response 
as a function of T or can be found in Fig. 1 for several hxed values of source 
M. In this case, when M = 0 there can be nontrivial solutions below a certain 
critical temperature Tc, which correspond to spontaneous symmetry breaking. For 
T > Tc and M = 0, the U{1)a symmetry is not broken and we have the same 
conhguration as [8]. When /is/T —)■ 0, the system becomes the previous case with 
/i = /is = 0. We have checked that the numerical solutions for the response in this 
limit satisfy precisely the relation (2.19). 



3 


10 


20 

/f5 

T 


30 


40 





Figure 1: The response O as a function of T and /is for different sources M. Left: Fixing 
T and B to be 8Ba/{7r‘^T‘^) = 0.1 and varying /is for qM/T = 0 (blue), 0.5 (purple), 1 
(red) and 2 (black). Right: Fixing /is and varying T and B at the same time with fixed 
ratio 8Ra/(7r^r^) = 0.1, qM/fi^ = 0 (blue), 0.04 (orange), 0.1 (red). 

• The third kind of simple while interesting solutions is for At = Vz = 0 while in 

general F), Az non zero (i.e. /is = 0, /i 7 ^ 0). In this case for M = 0 the solution for 

0 will be trivial despite /i 7 ^ 0 : no solution with spontaneous symmetry breaking 

will exist due to the zero axial chemical potential. The solutions of 0 for M 7 ^ 0 

are different from the solutions in the /i = /xs = 0 case due to nonzero Az for i? 7 ^ 0 

case and the solutions can be obtained numerically, which we will not give in detail 

2 

here. When 5 = 0 we have a simple solution V) = /r(l — (j-), I 4 = Aj = = 0 

and 0 is the same as (2.18). 

In the following, we will mainly consider the magnetoconductivity^^ for the first case. 
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Recently the magnetoconductivity has been widely studied in e.g. [25, 26, 27, 28, 29] for 
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i.e. /i = /is = 0 and will comment on the other two cases as well for Hall conductivity 
and the anomalous transport. A detailed study of the most general cases will be 
reported in a follow up work. 

2.2 Magnetoconductivity and relaxation time 

Without any axial charge dissipation, we know that the dual anomalous system will 
have an inhnite DC longitudinal conductivity in presence of a background magnetic 
held [8]. To show that our explicit breaking U{1)a model indeed encodes the axial 
charge dissipation, we need to compute the electric conductivity and the axial charge 
relaxation time. We consider the following huctuations at zero momentum on top of 
the background (2.12) 

, 6A^ = a^(r)e"*^*, <5$ = <hi(r, t) + i<h 2 (r, t) = (0i + . 

For the purpose of calculating the electric conductivity, we will impose the sourceless 
boundary condition for and 5A^ while 5Ez = iwvz at the boundary. We can always 
choose the gauge = 0 (or Ur = 02 = 0) using the fact that the equations in the 

bulk are invariant under the transformation 

6Vfj_ —>■ 6Vn + , 6Af^ —>■ 6A^ + d^A , <Fi —)■ <Fi — qA^2 ; $2 $2 + qA{(j) + $ 1 ). 

( 2 . 20 ) 

Note that in the gauge transformations above, the order of the gauge transformations 
should be at the level of perturbations, thus the terms qA ^2 and can be ignored 

as they are second order in perturbations. For these fluctuations at zero momentum, 
the equations for the longitudinal fluctuations a^, n*, 0i, 02 and transverse fluc¬ 

tuations ax,ay,Vx,Vy decouple from each other. In the /i = /is = 0 background, the 
fluctuation 0i also decouples from other longitudinal modes, which means that only 
the fluctuation of the phase mode 02/(g0) plays a role here. These equations can be 
found in appendix A.3.1 and A.3.2. We will focus on the longitudinal fluctuations to 
study the effect of axial charge dissipations on the longitudinal magnetoconductivity 
in this paper. The discussion on transverse magnetoresistivity and Hall conductivity 
can be found in appendix A.3.2. 

AdS 4 /CFT 3 with momentum dissipations. The longitudinal magnetoconductivity is generally mono- 
tonically non-increasing as a function of the magnetic field [30], i.e. positive magnetoresistivity. 
Though there can be exceptions, e.g. systems with paramagnetic impurities can have isotropic pos¬ 
itive magnetoconductivity, which is caused by B dependent scatterings. Here we wish to emphasise 
that the negative (anisotropic) longitudinal magneoresistivity is driven by axial anomaly and can be 
the smoking gun for the existence of axial anomaly. 
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For the case /x = /is = 0 (i.e. At = = Vt = Vz = Q), we have three independent 

equations of motion for a*, 02 as follows 


2^2 


SBa . 2g002 

wa' + ^^^Vz + 2iqr‘^f[ - 020 ' + <^(^ 2 ) 


. ^ 3 , 2g> 


„ , ,'3 , r\ , , a;2 ^ ^ 

* 0 + yJ 


r^f 


0 , 

0 , 

0 . 


( 2 . 21 ) 

( 2 . 22 ) 

(2.23) 


Note that the second equation is the equation of motion for and is of hrst 
derivative in the helds. The second order equation for 02 (i.e. the last equation in 
appendix A.3.1 with At = Az = 0) consistently follows from ( 2 . 21 ), ( 2 . 22 ) and the 
background equation of 0. One immediate observation is that S 5 = 0 which is dehned 
as the response of Jaz to the electric held Ez and this result is the same as the case 
without axial charge dissipation [ 8 ]. 

When r —)■ cx), we have 


( 1 ) ( 2 ) 
( 0 ) , 1 , 


( 0 ) I 1 I 

at = ay ^ - -inr H-- . 


(2.24) 


From (2.22), we have 


52 = W Co - 


iujir'^a't + SBaVz) 

2g2^5y02 


(2.25) 


where cq is an integration constant which can be chosen to be zero considering the 
sourceless boundary condition for 02. Thus when r —)■ cx), we have 


02 


g 0 



Inr 

2g2]Vf2 j-2 


+ 


iuj(2a[^'^ — %Bav^z^] 
Ar'^q^M'^ 


(2.26) 


From (2.21) and (2.23), we have — iojq^M ‘^cq and which 

show that the boundary constants a!'t\ ^ire fully determined by a® and Vt^\ Note 
that we always use the combination 02 /(q' 0 ) which is in fact the phase mode fluctuation 
—9{r) if we write 50 = 0(r)e*®*^^0 

To get the magnetoconductivity by solving (2.21 - 2.23), we impose infalling bound¬ 
ary conditions for Vz,at, 02 at the horizon and sourceless condition for at, 02 at the con¬ 
formal boundary, i.e. cq = = 0. At the horizon, the infalling boundary conditions 
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give 


V. = r 


’^o)" 


^(0) 


8-Baa(i)(4ro — ioj) + i;(o)Ci;(a; — 2iro) 


8rQ(2ro — ioj) 

at = {r - ro)“^ ( a(i)(r - tq) + • • • ) , 


(r - ro) + ... 


(2.27) 


/ 7rr[l/4]^(a(i)rQ(4ro — iu) + 32i?au(o)) 


^ 2 /(#) = (r - r„)--. ^ licfrlT[3H]HMITY 

where the dots denote terms which are higher order corrections. 


+ ... 


In numerics, notice that we have only two linearly independent solutions for (2.21 
- 2.23) which are fully determined by the near horizon values a(i),n(o). Using the fact 
that the system is invariant under the residual symmetry at —)• a^ + icjA, 02 —t 02 “ 5 'A 0 , 
where A is a constant independent of r, we will be able to generate solutions with cq = 0 
for each independent numerical solution. Then we can use their linear combination to 
set a® = 0 in numerical calculations [31]. 

We will identify the coefficient in front of 1/r^ of the fall-offs of Vz as the dual 

( 2 ) 

response, i.e. we are using the covariant current with {JzJz)r = We plot the AC 

longitudinal magnetoconductivity in Fig. 2. We hnd that the Drude peak behaviour 
becomes less obvious if we increase strength of U{1)a breaking M or decrease the 
magnetic field B. The appearance of the Drude peak is caused by the quasinormal 
mode with purely imaginary frequency u = —i/r^ approaching the real axis. Since we 
expect no new quasinormal modes to appear or vanish upon varying the magnetic field 
a sum rule of the form 

d r°° 

—— / Re[cr(a;, i?)](ia; = 0 , (2.28) 

d-B Jq 

is suggested to hold. By considering differences between AC conductivies at different 
values of the magnetic field we have checked numerically that this sum rule indeed 
holds. 


2.2.1 DC conductivity: negative magnetoresistivity 


From the hydrodynamic calculations in [8] the formula for the longitudinal magneto¬ 
conductivity in the = 0 limit is^^ 

e+p e+p 


cr = a E 


i {8Baf 

X5 


(2.29) 


^^Note that we chose the renormalisation scheme (2.17) such that the response is totally determined 
by the coefficient of the subleading term. 

^^In order that we can treat as an approximately “conserved” quantity, we should have T5T ^ 1. 
We also use the fact that the anomaly constant c = 8a for the covariant current. 
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Figure 2: The behaviour of both real and imaginary parts of AC longitudinal magneto¬ 
conductivity as a function oiu/T for different values of the source qM/T and the magnetic 
field SBa/n'^T'^. Top plots: 8B‘^= 0.1 with qM/T = 2 (black), 0.5 (red), 0.1 
(blue). Bottom plots: qM/T = 1, with 8B^a;^/(7r^T^) = 1 (blue), 0.5 (red), 0.1 (black). 
We found that the Drude behaviour becomes less obvious (or equivalently rs decreases) 
when we increase strength oiU{\)A breaking M or decrease the magnetic field B. Moreover 
we have checked numerically that the sum rule ^ Re[(T(a;, B)\dijj = 0 holds. 


where xs = denotes the static axial susceptibility and rs is the axial relaxation 
time. In the probe limit, as long as we have axial charge dissipation the resulting DC 
conductivity will be finite. The most straightforward way to see the effect of axial 
charge dissipation in the longitudinal magnetoconductivity is to go to the DC limit 
directly and in this limit we have 


{SBafr^ 

a — aE ^ - 

As 


(2.30) 


Note that this hydrodynamic formula only applies in the hydrodynamic limit t^T S> 1, 
however, holographic calculations below will apply even beyond this hydrodynamic 

limit. 


Instead of obtaining the DC conductivity by studying the AC conductivity and 
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then considering its cu —)■ 0 limit, there is an analytic way to directly obtain the DC 
conductivity following [32]^® using radially conserved quantities [33]. The idea is to 
perform the holographic calculations for conductivity in the exact DC limit, i.e. ca = 0. 
Then we can consider the following perturbations around (2.12)d® 


0, 0, -Et+V;,{r), 0), 0, 0, a^{r), a^(r)), Scj) = 0i(r)+i02 (r). 

(2.31) 

Substituting these perturbations into the equations (2.3 - 2.5) we obtain their equations 
which can be found in the appendix A.3.3. Note that these equations do not depend on 
the background helds Id, K- For zero axial density solution ps = 0, i.e. = 0, 

we have 


// , 3 , 2gV2 ^ ^ ^ ^ 

- 0 ) 


< + 1 - + y P.+ 




Near conformal boundary r —)■ cx), we have 


^( 1 ) ^( 2 ) „( 2 ) 
at = + ^lnr + ^ + ..., v^ = + ^ + 






(2.32) 

(2.33) 


(2.34) 


with a^t^ = From the equation (2.33) we have the conserved quantity along 

the radial direction J = —r^fv'^ — SBaat with drJ = 0, which means J|t —,-00 = -f|r-,-ro- 
When r —)■ CX), we impose the sourceless boundary condition for at, i.e. = 0. 
Thus the electric current which is the response of the external electric held can now be 
calculated at the horizon to be j = J|r^oo = J\r^ro- 

The next thing that we need to do is to get J\r^ro from the near horizon boundary 
conditions. At the horizon we have 


Vz = 

at 


E 


ln(r - ro) H-, 


2E{SBa) 

(47rT)g2r20o^ 




(2.35) 

(2.36) 


Note that near horizon the coefficient of the subleading term in at is a free parameter 
which is precisely the shooting parameter that can be used to determine the sourceless 
condition for at near boundary. Also note that the near boundary condition (2.36) 

^^For a recent application and development of this method see e.g. [25, 27]. As will be shown in 
the subsection 4.2, this DC method in our case is equivalent to directly taking the a; —)■ 0 limit for 
small oj AC conductivity obtained from near-far matching calculations. 

^®Note that one can also work in the radial gauge Or = 0 and the same results should be obtained. 
Here we include in order to totally ignore equations related to (f> 2 - 
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seems to be inconsistent with (2.27), however, this is dne to the different gange: if we 
work in the gange (/)2 = 0, we shonld keep and the leading order in (2.27) shonld be 
the same as here in the limit ca —>■ 0. It follows immediately 


ctdc = = vtT + 


E ' 7r3T3g2(^2 


(2.37) 


with 00 the horizon valne of the backgronnd scalar field. We shonld emphasise that the 
above expression does not depend on the mass of the scalar filed in the bulk and it also 
applies for the case /x 7 ^ 0 as long as ps = 0. We emphasise that the dependance 
is exact, and the recent experiments [14, 15, 16, 17] also found this B^ dependence in 
the longitudinal magnetoconductivity. 


For the simple case of /i = 0, i.e. Vt = Vz = At = = 0, (p has an analytic solution 

(2.18) for mj = —3, and we have (2.19), thus 


o'dc 


_ r[i/4p (sag/ry , 

T ^ 4ir2r[3/4]2 (qM/TY ^ ^ {qM/TY ' 

This is one main results of our paper. In fact this DC magnetoconductivity can also 
be obtained analytically using a near far matching calculation for small u, in which 
situation the near horizon boundary conditions (2.35) come from solutions in the near 
region with infalling boundary conditions while the equations in the DC limit are 
exactly the far region equations. In the subsection 4.2 we show the near far matching 
calculation which produces the same DC result for both this explicit breaking case 
and the massive gauge field case. From the near far matching analysis it is natural 
to identify ttT as the quantum critical conductivity aE in (2.29) and (2.30) in this 
^ -C limit, and the value of aE strongly depends on the value of We will 

explain at the end of the next subsection and also in subsection 4.2 that the discrepancy 
between this ag = ttT and the ag in the case without axial charge dissipations (A.2) 
in [ 8 ] comes from the noncommutative nature of the two limits ca —)■ 0 and Ts — )■ cxd. 


From the formula (2.38) for the DC longitudinal magnetoconductivity, we can see 
that the negative magnetoresistivity has a universal B^ behaviour even in the large B 
limit (quantum regime). In the regime where the axial symmetry breaking is small this 
can be understood from the B dependent behavior of the charge relaxation time rs such 
that Ts/xs in (2.30) does not depend on B which results in the universal B^ behavior 
in (2.37). Our numerical results obtained from the AC conductivity with ca —)■ 0 also 
match quite well with the formulae (2.38) as can be seen in Fig. 3. 


2.2.2 Calculating Xb rs independently 

With the analytic result (2.38) for the DC magnetoconductivity, we cannot identify 
Ts and Xb because there is not enough information in this formula. Thus we will 
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Figure 3: The dotted lines are the numerical results for the DC conductivity obtained 
numerically (by setting uj/T = 10“® in the AC conductivity calculation). Note that if t^T 
is large, we should choose extremely small oj/T to obtain the DC conductivity numerically. 
The solid lines are the best fitting for numerical data. Left plot is for qM/T = 1 and 
the best ht is a/T — vr = 0.00005 + 21.600(8i?Q:/(7r^T^))^/(gM/r)^. Right plot is for 
8Ba/{TT^T^) = l/vr^ and the best fit is a/T-ir = 0m07+21.574{8Ba/{TT^T^)f/{qM/Tf. 
The fitting formulae exactly reproduces (2.38) from analytical calculations. 


calculate the static susceptibility X 5 and the charge relaxation time Ts independently 
and substitute them into (2.30) to be compared with the analytic result (2.38). 

The static susceptibility xs = limtj_^o{'^A'^A) ^an be calculated from perturbations 
(2.31) while setting E = 0, and the boundary conditions are now at,Vz being regular 
at the horizon and being sourceless at the boundary, the latter of which can be fixed 
from gauge transformations. We need to solve the following equation 



and this equation can only be solved numerically. The result is shown in Fig. 4. Note 
that different from the DC conductivity which can be totally determined by the near 
horizon data, the static susceptibility cannot be fully traced to horizon data which is 
due to the fact that U{1)a is no longer an exact symmetry in this explicit breaking 
model. 

The next step is to calculate independently. This can be done by analyzing the 
quasinormal modes at fc = 0 on top of the /x = /is = 0 background with no source 
for the electric held. We will solve the {vz,at,4>2) in (2.21 - 2.23) in the complex 
frequency plane with infalling near horizon boundary conditions for atiVz-i(l) 2 '- (2.27) 
with Co restored in 02 - As stated above there will be two linearly independent solutions 
for Co = 0 with boundary values qMc^) and qMc^^), which are 
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Figure 4: Left: X5 a-s a function of B for different sources with qM/T = 0.1 (black), 1 
(red), 2 (blue). Right: X5 as a function of M for 8Ba/{TT‘^T^) = 5 (black), 2 (red), 0.1 
(blue). 


oj dependent. The residual gauge symmetry will give the third solution of boundary 
value (0, —iojA, qMA) with radially independent A, and this represents the degree of 
freedom due to Cq. We can define the matrix for boundary values as 




af' 

qMc^ ^ 

Mbnd — 



qMc^^ 


\ 0 

—iujA 

qMA j 


Following [31], the QNM frequency is given by the zeros of the determinant of the heled 
values at the boundary, i.e. |Mbnd(—^a;/)| = 0. We define the axial charge dissipation 
time Ts as Ts = l/uj. This is consistent with the hydrodynamic modes following the 
hydrodynamic equation = 0 and = —^Ja [ 20 ] although for small t^T 

these equations may not apply anymore given that other non-hydrodynamic QNMs 
will dominate. 

The hnal numerical plot on is shown in Fig. 5. One can see that t^T increases 
when we increase the magnetic field or decrease the value of the source that explicitly 
breaks U{1)a- When t^T 3> 1 hydrodynamics applies^^ and we are expecting a Drude 
peak behaviour [ 8 ] (i.e. coherent metal behaviour) and this is exactly what we found in 
the small frequency regime of the AC conductivity. When t^T < 1 , the hydrodynamic 
description breaks down (incoherent metal behaviour) and the contribution from this 
hydrodynamic QNM will not be dominant [35]. Consequently there is no Drude peak 
in this regime in the AC conductivity. Moreover the standard Boltzmann theory which 
is based on the the quasiparticle picture does not apply for the small case. Here we 

this case in principle could be determined by memory matrix formalism as the three dimen¬ 
sional case [34]. We leave this interesting question for future investigation. 
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want to emphasize that even when the axial symmetry is strongly broken we still have 
negative magnetoresistivity as we have shown in the previous subsection. We note that 
recently in a weak coupling context in [36] it was pointed out that the ionic scattering 
can produce universal positive magnetoconductivity (or negative magnetoresistivity) 
with behaviour for a generic 3D metal in presence of parallel electric and magnetic 
fields. 



Figure 5: Left: t^T as a function of B for different sources with qM/T = 0.1 (black), 1 
(red), 2 (blue). Right: t^T as a function of M for SBa/in'^T'^) = 5 (black), 2 (red), 0.1 
(blue). 


2.2.3 Scaling behaviours of X 5 and rs on B and M 

We will analyse the scaling behaviour of X 5 and with B and M in this subsection 
and we will show that in the hydrodynamic limit the hydrodynamic formula reproduces 
the dependence of the DC conductivity on B. First, note that the formula 

o^DC = o'E -\ — -{SBa)^ (2.41) 

X5 

for the DC conductivity is only applicable in the hydrodynamic limit B/T‘^ 1 and 

t^T 1 which means M/T 1. However, in the following, we will show the scaling 
behaviour of X 5 aad rs on B for a large range of B, in which this formula still coincides 
with the holographic DC result. This also happen in the case without any axial charge 
dissipation [ 8 ]. 

Fig. 6 shows that when B is large (fixing M, T), both X 5 and rs are linear in 
B, while when M/T is small (fixing T, H), xs is a constant and ts (or F 5 = l/rs) 
is proportional to M~^ (or M^). Fig. 7 shows that when M/T is small, r^/Xb does 
not dependent on B in the large B regime. From these hgures we can see that the 
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formula (2.41) is valid even for large magnetic field in the regime of small for small 
M/T. When t^T S> 1, we have erg = ttT and Ts/xs = 21 .6/(7r^T^g^M^). At small M, 
'I'b/Xb is proportional to M~^. 




M M 

T T 


Figure 6: Scaling behaviour of static susceptibility X5 and relaxation time Here the 
prime means the derivative of the parameter in the horizontal axis. Black solid line is for 
M/T = 0.1 while blue solid line is for = 0.1. 


Now we explain why as that we observe here is different from the one in the case 
without any axial charge dissipation in [8]. Note that if we take the limit rs —)• cx) 
(by M/T —)■ 0) in the DC conductivity (2.37), we can see that the quantum critical 
conductivity aE = ttT in this limit and we cannot see the effect from the anomaly 
and the magnetic held as in the case without axial charge dissipation [8], where 
drops with B at small B and approaches zero at large B (see (A.2)). This seeming 
discrepancy comes from the non-commutative nature of the two limits of a; —)■ 0 and 
Ts —)■ oo. In the case without axial charge dissipation, the limit of Ts —)■ cxd was taken 
hrst while in our formula for the DC conductivity with axial charge dissipations, the 
limit of a; —?• 0 was taken hrst and some terms in a e disappear in this limit. This 
means that cte in the formula in [8] can only be seen in the limit l/rs <C cn — )■ 0 .^® 

^®A similar phenomenon was observed in the momentum dissipation case. For four dimensional 
Reissner Nordstrom black hole solution without momentum dissipation, as behaves nontrivial at low 
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Figure 7: Left: t^Ix^ as a function of B for different sources with qM/T = 0.1 (black), 1 
(red), 2 (blue). Right: rs/xs as a function of M for SBa /= 5 (black), 2 (red), 0.1 
(blue). Orange dashed line is the constant coefficient 21.6 in Eq. (2.38). This is consistent 
with the fact that our formula Eq. (2.41) is only valid when T 5 T ^ 1. 


2.3 The anomalous transport coefficients 


Another interesting transport property for the explicit axial symmetry breaking case 
is that we can still observe a similar anomalous transport as the chiral magnetic effect 
as in the normal chiral anomalous systems. To show the anomalous transport in this 
case, we should consider perturbations of the system at a; = 0 and small k and use 
the Kubo formula as in normal chiral anomalous systems. We will consider the chiral 
magnetic conductivity J = cr^B, the chiral separation conductivity J5 = ctcseB and 
the axial magnetic conductivity J5 = assBs, where B5 is an axial magnetic held. 

We consider the transverse huctuations ax,ay,Vx,Vy with k along the z direction 
and their equations can be found in the appendix A.3.4. It is convenient to introduce 
a± = Ox ± iay,v± = Uj, ± iVy. For the /i = ps = 0 case 47 = W = A = ^2 = 0, it is 
easy to hnd that as = <755 = ctcse = 0 . 

For the case with nontrivial axial charge density (i.e. /i = 0, /xs 7 ^ 0) we have 


J' 


r 




, 3 f 

■± +1 ^ j 

,2 


+ 




(2.42) 


Since the vector sector and axial sector decouple from each other, ucse = 0. 

Up to order A:, the asymptotic series at conformal inhnity is = Uq + ^ + • • • cind 

temperature which can be found in [37]. However, after introducing momentum dissipation through 
different mechanisms [38, 39, 25] as is constant. The reason should be the same as our case: it is 
exactly due to the non commutativity between two different limits Tm ^ 00 and a; —>■ 0 . 
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(2.43) 


a± = ao + ^ Inr + ^ + ... with ai = —a^q^X^. The solution for v± is 


v± 



Note that T^G^y + = G± and as = hnifc^o = hmfc^o ‘^^ 2 k ■ = 

Sa/is- We note that the exactly conserved (consistent) current that is dehned via the 
fnnctional variation of the on-shell action with respect to the gange potential has an 
additional contribution from the Chern-Simons part of the action = Aae^’^^^A^.J^px. 
This part is easy to calculate since it is completely determined by the bonndary valnes of 
the helds and adds —Saps to the chiral magnetic condnctivity. Therefore we hnd as = 
Sa/is for the current without the contribution of the Chern-Simons term {covariant 
current) and o'B(con) = 0 for the exactly conserved current. This is in line with the 
expectations from recent arguments [ 6 ] that exactly conserved currents cannot have a 
non-vanishing expectation value in eqnilibrium.^® 

Qa _ Qa 

Similarly (J 55 = lim^^o ^ 2 fc with the two point correlator for axial gange held 
a±. We can expand a± = -|- ko!'^ -|- + ... with 


( 0 )" /3 /' 


( 0 ) 

af 


= 0 , 


8 a 


= T 


r^f 




(2.44) 


Since one can rescale a||!^ —)■ apsa^^ we conclnde that (T 55 is proportional to a/xs. 
The exact result can only be obtained nnmerically, which is shown in the left plot of 
Fig. 8 . We note that in the limit of large mass M//r 5 , (J 55 approaches |a/i 5 . Also 
note that this behaviour should not depend on B for hxing SBa/Ti'^T'^. This is the 
same nniversal valne that was observed in anomalous holographic superconductors in 
[21, 22, 40] in the T —>■ 0 limit. We also note that for low temperatnres we end np 
in the M —)■ 0 limit in the snperconducting phase. Therefore the bine line in hgure 
8 varies very little. In order to interpret this resnlt we again add the Chern-Simons 
contribntion to the current in order to obtain the (consistent) cnrrent that conples 
axial gange held Ay. This time the Chern-Simons current is 

The factor 1/3 can be nnderstood by noting that this comes from the triangle anomaly 
with three eqnal current on the vertices {J 5 J 5 J 5 ) which implies a symmetry factor of 
1/3 compared to the triangle with one axial and two vector like currents. The Chern- 
Simons term contributes now — This means that in the large mass limit M —)■ cx) 

^®For a detailed discussion of covariant vs. consistent definition of currents in relation to anomalous 
transport see [5]. 
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in which the axial symmetry is maximally broken the total axial current va ni shes! This 
seems a very intuitive result. 

Finally we consider the case /xs = 0,/r 7 ^ 0. For simplicity we will discuss the 

chiral separation conductivity in the linear response approximation, i.e. with vanishing 

background magnetic field. When B = 0, we have a simple solution with non vanishing 
2 

fields Vt = /r(l —and 0 as the same as (2.18). In this case as = <755 = 0 while nonzero 
fcsE which can be found in the right plot of Fig. 8 . Let us make a comment on the 
behaviour of the chiral separation conductivity. In this case there is no contribution due 
to the Chern-Simons term to the current since we only switch on a chemical potential 
for the conserved vector like symmetry /x 7 ^ 0. We find that also in this case the axial 
current induced by a magnetic field vanishes in the limit of maximal axial symmetry 
breaking. Since in this case there is no Chern-Simons current also the covariant current 
vanishes. 



qM M 

T 


Figure 8 : Left plot: <755 as a function of the source M for SSa/vr^T^ = 0.1, a = 1, 
F/ps = 0.05 (blue), 0.075 (red), 0.1 (black). Note that Tc/ps ~ 0.765 for this case. 
Dashed line £ 755 /( 0 ^ 5 ) = 8/3. Right plot: ucsE as a function of M/T for R = 0 and 
P 5 = 0 while /X 7 ^ 0. When M —>■ 0, we have ctcse = 8 ap. For large M, cxcsE —>■ 0. 


The important conclusion of this analysis is that in the limit of maximal axial 
symmetry breaking via the mass parameter M the expectation value of the axial current 
J 5 vanishes for both the chiral separation effect and the axial magnetic effect, but only 
if one uses the consistent definition of the currents. Since a vanishing axial current for 
maximal axial symmetry breaking seems a physically plausible result we take this as 
an argument in support of using the consistent definition of currents. 
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3 Holographic massive U{1)a x ^(1)v model 

In this section we will concentrate on the massive U{1)a model and consider the pres¬ 
ence of two U{1) gauge fields in the bulk. Both fields are coupled via a Chern-Simons 
term as in the previous model. In order to implement (axial) charge dissipation we 
consider a constant mass term for one of the U{1), which we refer to as the axial U{1)a- 
This is achieved in a gauge invariant manner via the Stiickelberg mechanism [19, 20]. 
The action reads 

5 = I (fi + tl) - ^ (.4, - s,e) (/I- - a'-e) 

- a„«) (F„gF^s + SF^gFgi) ]) , (3.1) 

where F = dA and F = dV. The Stiickelberg field and the gauge field transform 
as 9 —)■ 6* -|- A, —)■ leaving the mass term invariant. The Stiickelberg field 

it is not charged under the U{l)v- 

As in the previous section we will work in the probe limit with Schwartzschild AdS 
background. The equations of motion are 


- m\A'^ - d^e) + = 0 , (3.2) 

+ 2ae^^^^PF^0Fyp = 0 , (3.3) 

Vp {A>^ - d>^e) = 0 . (3.4) 

The mass term affects the asymptotic behaviour of the axial gauge field A^. The 

asymptotic expansion of the different fields reads 

ApAp{Q)r^ + ... + Ap(Q)r ^ ^, (3-5) 

— k)i(o) + ••• + ^Ai(o)^ ^ , (3.6) 

0 ~ 0(0) ...0(o)r“'‘ -I- ..., (3.7) 


with A the conformal dimension of the source for dual axial current operator A = 
—1 -|- \/im? F 1. A^(o), V)j(o), 0 ( 0 ) are the coefficients of the non-normalisable modes and 
^( 0 ) are the coefficients of the normalisable modes. As one can see the 
presence of the mass changes the dimension of the operator dual to the axial U{1)a'- 
[J^] = 3 -|- A. This implies that the dual axial charge is not conserved since in a 
conformal theory a conserved current must saturate the unitarity bound. We refer to 
A as the anomalous dimension of the axial current. Moreover the anomalous dimension 
has a bound, given by the condition A < 1. This is obtained by requiring the dual 
axial operator to be irrelevant in the UV. 
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It is necessary to renormalise the theory in order to obtain hnite observables [41], 
The amount of divergent counterterms needed depends on the precise value of the 
mass and diverges as one reaches the marginal case A = 1. For this reason we will 
only consider masses such that A < 1/3, which keep this amount minimal. For these 
values of the mass the boundary term containing the counter terms reads 








(3.8) 


with A comment is in order regarding the Ward identities of the 

currents. The vector current is conserved and therefore the divergence vanishes 


= 0 . (3.9) 

With such a choice one could expect the divergence of the axial current to be explicitly 
proportional to a A J- + F A F). However this intuition fails in the massive case. 
The axial symmetry is broken by the dynamical internal SU{N) degrees of freedom. 
This implies that the “current” operator is to be considered as a non-conserved current 
which, therefore, lacks any constraint given by the symmetry. This is nicely seen in 
this model by computing the expectation for the divergence of this current. With the 
renormalised action and after using the asymptotic expansion one gets 

d^J>^ = 0, = (2 + 2A)a^i5,), (3.10) 

with no explicit constraint for the axial current. 


3.1 Magnetoconductivity and relaxation time 

In this section we compute the electric DC conductivity, the static susceptibility and 
the axial charge diffusion time in presence of a background B held. As background 
we switch on a spatial component of the vector held 14 = By which trivially fulhls 
the background equations of motion and generates a constant magnetic held aligned 
to the 2 ; direction. To compute the appropriated quantities we switch on perturbations 
on top of this background with hnite frequency and momentum aligned to the B held 
59 = 5A^j_ = and 51// = For our purposes 

we can just focus in the sector that contains 02 , , n*, r/ which is decoupled from 

the other components in our background. The explicit form of the equations can be 

addition to the divergent terms we have implicitly chosen a scheme that respects the axial 
symmetry by including the hnite term in (3.1). 
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found in appendix A.4. Here we list the fc = 0 sector, 


+ ^3 


, lujm 

V. H- —V = 0 , 


rV 


,, 3\ , a;2 SaB , ^ 

{j + -jvz + ’ 

9 9 /> / / / \ 

—im r fif] + uj[a^-\ -= 0 . 


(3.11) 

(3.12) 

(3.13) 


In order to obtain two-point functions and quasinormal modes numerically^^ we 
build the numerical, matrix valued, bulk to boundary propagator F with the appro¬ 
priated normalisation such that 


/r-^ai{r)\ / am\ 

v,{r) =F(r) u,(o) , F(A)=I, (3.14) 

\ V{r) / \h(o)/ 

where A is the cutoff radius. This can be obtained imposing infalling boundary condi¬ 
tions and a set of orthonormal values of the amplitudes at the horizon. With the usual 
holographic prescription different correlators can be written as a linear combination 
of derivatives of F(r) and quasinormal modes can be obtained from the zeros of the 
determinant of F~^(r). Concretely for the electric DC conductivity and the axial static 
susceptibility we hnd 


= lim A3+^F1^,,^(A) 

u=k=0 

ctdc = lim lim A^— (F^ ^^^(A) + a;^F^^,^^(A) log(A)) 

A->.ooaj-s>0 lOJ 

where subscripts at, Vz refer to the appropriated entry of the matrix. 


(3.15) 

(3.16) 


3.1.1 Electric DC conductivity 

We compute the negative magnetoresistivity in this model and compare it to the results 
in the previous model. Note that the anomaly term proportional to the external electric 
and magnetic helds cannot be implemented in the hydrodynamic expansion in this 
model. This can easily seen from the conservation law since Jg-|-cE-B , does 

no longer hold due to the anomalous dimension of the axial current. This ultimately 
implies that one cannot derive an equation analogous to (2.29) in this situation. 

refer the reader to [20] for a thorough explanation on how to numerically compute different 
quantities within this model. 
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Figure 9: Dots show numerical result for the DC conductivity obtained by setting 
oj/T = 10“^ for the AC conductivity. Solid lines show the best fit for the data. Left; 
Fixed m = 0.5 and the best fit is a/T — vr = 0.0001 + 0.99947r(8ila/7r^T^)^/m^. Right: 
Fixed SBaKir'^T"^) = 4 and the best fit is a/T — tt = 0.0001 + 0.99947r(8RQ;/7r^T^)^/m^. 
The fitting formulae exactly reproduces (3.22) from analytical calculations. 


As we will show now the magnetic held dependence of the DC conductivity in this 
model still implies positive magnetoconductivity. We hrst obtain analytic formula for 
DC conductivity by means of the near horizon analysis [32]. The procedure is the same 
as the one in previous subsection 2.2.1. Consider huctuations 


SV^ = {vt{r),0,0,-Et+v^{r),0), 6Af, = {at{r),0,0,a^{r),ar{r)), 69 = r], (3.17) 

we need to consider 


// 3 , m" 

H—m-+ 


8aB 


r^f 


-V, = 


r 




+ 


8aB 

r^f 


■m = 


0 , 

0 . 


(3.18) 

(3.19) 


Near conformal boundary r —)■ cx), we have at = +... ,Vz = + 

+ .... From the equation (3.18) we have conserved quantity J = —r^fv/—8Baat 
with drJ = 0. Thus Jlr-s-oo = J\r^ro- When r —)■ CX), we impose the sourceless boundary 
condition for at, i.e. = 0. Thus the electric current which is the response of external 
electric held j = J|,._).oo- We are looking for the solution which is regular near horizon. 
Thus we have 


Vz 


E 

— ln(r - ro) + 

4E(,SBa) ^ 

(47rT)m^rQ 


r - ro), 

ro) ■ 


(3.20) 

(3.21) 
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Note that near horizon the subleading term in at is a free parameter which is precisely 
the shooting parameter that can be used to determine the sourceless condition for at 
near boundary. It follows immediately 



(3.22) 


This formula shows again the key point of this work: the DC longitudinal magneto¬ 
conductivity depends quadratically on the magnetic field independently of its strength 
when charge relaxation is built in the models. It is remarkable that this behaviour holds 
in this setup too, despite to the lack of a clear hydrodynamic prediction in terms of 
and X 5 - We also checked the formula numerically. Our results show perfect agreement 
with the analytic formula, as shown in Fig. 9. 

The AC conductivity in this model has been studied numerically already in [20]. 
Here we note that also in that model a sum rule of the form (2.28) was found to hold. 

3.1.2 Axial charge dissipation time 

We compute the dissipation time for the axial charge from the gap in the imaginary 
part of the lowest QNM associated to this symmetry. As shown in [20] and in [42] in 
chiral kinetic theory a diffusive mode for the vector charge and a gapped dissipative 
mode for the axial charge. It is this gap that we identify with the inverse of the axial 
relaxation time. 

We compute the QNMs at zero momentum numerically using the numerical tech¬ 
niques explained in [31]. The QNMs are obtained from the zeros in the determinant 
of the inverse bulk to boundary propagator defined in the previous subsection (3.14). 
The explicit equations are shown in appendix A.4. In hgure 10 we show our results 
for the relaxation time. As shown in [20] the relaxation time is found to be inversely 
proportional to the bulk photon mass. Moreover the behaviour of this relaxation with 
increasing magnetic held shows a transition to a linear regime for big enough magnetic 
helds, see Fig. 10 for £t analysis. This is analogous to what was found in the explicit 
U{1)a breaking model. 


3.1.3 Static axial susceptibility 

Although (2.29) does not straightforwardly apply to this model too it is still interesting 
to compute the axial susceptibility and its dependence with the magnetic held. 
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Figure 10: Left panel: Axial charge relaxation time as a function of magnetic field for 
different values of the bulk mass m = 0.15 (Black), m = 0.25 (Red), m = 0.5 (Blue). Right 
panel: Numerical data (dots) and linear fit (red) for at m = 0.5 and large magnetic 
field. The best fit in this region is t^T = 0.0979 + 0.5030^^^^^. 


To obtain the analytic result we set E = 0, then we have 


„ 3 , 


at = 0. 


(3.23) 


Near horizon we impose at{ro) = 0. The solution is 

^ _ ^2/2] 


at = i[ — 


T [1 - Vl + my4 + ^1 - /3V4] T [1 + + ^1 -/3V4] 

Fi [(-a/1 + + ^1 - /32)/4, {Vl + - v^T^^^)/4,1 - \/l - /2,r'^/r^] 

_ ^2/2] 


X 




r [1 + Vl + m2/4 - v^i - /^V4] r [1 - Vl + my a - ^i - /3V4] 

2 F 1 [(-Vl + m2 + a/ 1 - /92)/4, (Vl + m2 + ^1 - /92)/4 ,1 + V/r^] , 

where /3 = Near conformal boundary the above solution behaves as at = + 


S ) —1—vT-j-t 


thus 


A5 


= -2VI Tm^V 


4 -) 


(+) 


= 4r, 


r[i-^]r[i(^r+^-yr^)]r[i(vr+^+v/r^)] 


Now let us study the behaviour of static axial charge susceptibility in large B limit. 
When f5 = —>■ 00 , we have 


A5 ^ 4ro 


2vTTr^ 
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\/l+T1 
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(3.24) 


































































Figure 11: Left: Static susceptibility as a function of magnetic field for different values of 
the bulk mass m = 0.15 (Black), m = 0.25 (Red),m = 0.5 (Blue). Right: X5 for ru = 0.5. 
Red line corresponds to linear fit and blue line to {8aB /j'2^2^g^g expected from 
(3.24). 


We check this numerically. We compute the axial static susceptibility X 5 and its 
dependence with the background magnetic held numerically by means of the the Kubo 
formula 


= (yy) 


(3.25) 


uj=k=0 


In hgure 11 we show the behaviour of the static susceptibility against magnetic held. 
As expected there is a transition to a hxed scaling for large magnetic held. However, 
as predicted in (3.24) the exponent is now V^lTT^^, see Fig. 11 for ht analysis. Our 
results show indeed that the ratio r^/Xb is not B independent for large values of 5, 
contrary to what was found in the previous model. 


4 Remarks: relations of two models and DC results 
from small uj matching 

In this section we will make some important remarks on two separate aspects: the rela¬ 
tions between our two holographic axial charge dissipation models and rederivation of 
the DC conductivity based on radially conserved quantities from the near-far matching 
calculation for AC conductivity at small frequencies. 
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4.1 Comparison between the two models 

These two axial charge dissipation models look very different: one breaks the charge 
conservation symmetry by giving a mass to the gauge held and the other breaks this 
symmetry explicitly by a scalar operator source. However, if we look closer to these 
two models we will hnd that these two are in fact closely related. We can compare the 
actions or the equations of motion of these two models. It is easy to hnd that in the 
equations of motion for perturbations (3.11 - 3.13), if we replace in the massive 
gauge held model by and rj ^ we will reproduce exactly the same equations 
of perturbations for the explicit breaking model (2.21 -2.23). It explains that our DC 
conductivities (2.37) and (3.22) in these two models have a universal formulae. This 
shows that the scalar operator in the explicit breaking model gives an ehective mass 
to the gauge held.^^ By solving the equations of the two models, we can even see that 
the massive gauge held model is in fact a special case for the explicit breaking case in 
which the bulk mass for the scalar is set to zero and p = /is = 0. The massless scalar is 
dual to a marginal operator. A source for it does however change the scaling dimension 
of the axial current. An advantage of the explicit breaking model is that we can choose 
the mass of the scalar held freely without changing the the scaling dimension of the 
axial current. 

Note that in both models we can choose two kinds of gauges (or any combination 
of the two): one is to choose SAr = 0 and the other is to choose 5r/ = 0 in the massive 
gauge held case or 64>2 = 0 in the explicit breaking case. For both these two gauges, 
we will hnd that the equations for perturbations are the same with the mass of the 
massive gauge held replaced by an ehective mass generated by the scalar source. 

Finally, though at the level of actions or equations of perturbations the two are 
closely related or equivalent in some sense, there is a subtlety here that in the massive 
gauge model, there is no Higgs mode, i.e. r/ is a real scalar while 0 is a complex 
scalar held. This Higgs mode does not have any ehect in transport coefficients, but it 
may have other ehects. In particular when /r 7 ^ 0 or /is 7 ^ 0 the Higgs mode will not 
decouple. 

4.2 Near far matching calculation for the AC conductivity at 
low frequency 

In this subsection we reproduce the DC magnetoconductivity result which was obtained 
using radially conserved quantities from the near far matching calculation for both the 

^^Similar physics happens for the momentum relaxation case: as concluded in [43] the explicit 
momentum breaking by a scalar field gives graviton a mass. 
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explicit U{1)a breaking and the massive gange field case. As we emphasied in the 
previons snbsection, the eqnations (2.21 - 2.23) and (3.11 - 3.13) are the same if we 
replace m? in the massive gange field model by 2q‘^(j)‘^ and h In the following 

we consider Eqs. (3.11 - 3.13). It is straightforward to apply to the replacement to 
(2.21 -2.23). 

We work in the gauge Srj = 0 (or 502 = 0 gauge for the explicit breaking case) 

r2 

and in the coordinate u = ^ for convenience. The equations for the perturbations 
6At = ai(M)e “*^0 514 = Vz{u)e~'^^^, and SAr = ar{u)e~^‘^^ are 

4m^(1 — M^)(4ai?n(. -t- rla” + ir^ua'^) — rn^rlat = 0 , (4.1) 

(—m^rg(l — u^) + rluu'^)ar — iuu{AaBvz + r^a^) = 0 , (4.2) 

4m(1 — u^) ^rQ((l — u^)v'^y + AaBa^j + uPvz + lQiaBuoj{l — u^)ar = 0 (4.3) 

with ' the derivative in u. 


The near region is defined as 1 — u <C 1 and the far region is defined as 1 — u 3> ^. 
We will first solve the near region and the near horizon boundary conditions for the 
far region are provided by the near region solutions expanded at the matching region. 

In the near region, we have the infalling boundary conditions and the solutions at 
leading order are 

n, ~ (l-n)-'‘"/""°(1+ ...), (4.4) 

at ~ (1 - (caso + si(l - a) -1- ...) , (4.5) 

a,. ~ (1 - { 0 JS 2 + 53(1 - m) -h ...) , (4.6) 

where sq, si, S 2 , S 3 are constants which may depend on uj and m and using the equations 
there will be only one free parameter which we denote as S 2 . The ’ above denotes 
subleading terms which contribute to the same higher orders in the equations. 


The following matching calculations are equivalent to defining new functions as 
Vz = (1 — at = (1 — (u) and a^ = (1 — and 

solving these new functions with boundary conditions at = 0 at m = 0. Here we take the 
terminology of matching in order to present the results more clearly to be understood. 
With the near region leading order solutions, we then expand them at leading order in 
uj in the matching region ^ 1 — u 1 , which can be used as near horizon boundary 

conditions for the far region: 

~ 1 - ^ln(l- m) , (4.7) 

4ro 


at ~ -4roS2Co’ -t- (1 - u) ( 


AaB 


2im^rls2) + + ^m^roS2)a; 


/ NT/ \ f 4o:7? Q 9 

- 1 - (1 — m) ln(l — u) ( —2 - 2im rQS 2 

V ^0 


—too 

4rn 


(4.8) 
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where S 2 is the tuning parameter to make sure that the boundary condition at = 0 at 
the boundary u = 0 is satished. is decoupled in the far region and is not important 
here. 

To derive the far region equations, we can drop all the terms in the equations 
at order o{u) while order 0{u) terms should be kept. In the equations above, we 
can solve from the second equation and substitute it to the hrst and the third 
equations and find that those 0(uj) order terms all become 0 ( 0 }“^) order terms and can 
be ignored in the far region. Note that in this procedure we have secretly assumed 
that m^rQ/a; S> uu/{l — u) in the far region as can be seen from the coefficient of 
in the second equation. With a further constraint that the terms in (4.1) and (4.3) 
are at order o{u}), we know that the following far region equations are only valid for 
S> cj/ro (this condition should be substituted by 3> oj/tq in explicit breaking 
case), i.e. cars <C 1: 

((l-«>;)' + ^a; = 0, (4.9) 

4n^(l — u^) ^— rn^at = 0 , (4.10) 

which are the same as we are studying in the paper (2.32 - 2.33) and (3.18 - 3.19) for 
DC conductivity and the hrst equation can be integrated to give 


(1 - + 


AaB 


O't — Cq 


(4.11) 


where Cq is an integration constant which can be decided by the near horizon analysis 
to be 


^0 = 


iu 

2ro 


16aB 

- S2UJ, 

ro 


(4.12) 


at order 0{u;). 


The next step is to see which S 2 can set the leading order coefficient of at at the 
conformal boundary to be 0. If S 2 has 1/uj dependence at the leading order, then there 
is only one term (1 — u) in at at the order l/u due to the term in the equation for 
at, which cannot be canceled by other terms at the boundary. Thus S 2 ~ S2o + i^S2i • ■ ■. 
Then for small uj we can see that there is one term (1 — u) in the near horizon expansion 
of at at order 0{1) while others are at 0{u). 0{u) and 0{1) terms belong to linearly 
independent solutions and at the boundary they will lead to boundary values at order 
0(u) and 0{1) separately because the far region equations do not depend on u. To 
make sure that the leading order in a* is 0 at the conformal boundary, we have to impose 
that the Oil) order coefficient in front of (1 — u) vanish, which gives S20 = — S21 

is now the tuning parameter to shoot the boundary value of a* to 0. 
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From the far region solution we know that v'^{u = 0) = Co and = 1 + ihiu/m?' 
where hi should be a constant at order 0{1) in both u and w?. Note that in ^^(0) there 
should be order corrections which can only be obtained by considering the subleading 
order equations in the far region. Thus the final result for the magnetoconductivity 
should be 


cr = 


2rQ Cq + h20j‘^ 


lU 


ro + 


0 


ih20J 


1 + ibi-^ 

i- rrt^ 


(4.13) 

(4.14) 


where 62 comes from the corrections to n^(0). With a simple estimate 62 ~ 1/m^ for 
small m. This is consistent with the fact that we are working in the small u:/m?T limit, 
in which the subleading order corrections compared to the leading order coefficients in 
the numerator and the denominator are usually at the same order and hi can not be 
identified as related to the relaxation time Note that the terms in the numerator 
at order 0{u^) is accurate and non-perturbative in First we can see that when 
a; —)■ 0 we reproduce the DC results: Eq. (2.37) in Sec. 2 for the U{1)a explicit 
breaking model and Eq. (3.22) in Sec. 3 for the massive gauge held case. Second, 
it is easy to see that when we take hrst a; —?• 0 and then 0 limit, ro will give 

the value of cte. In fact there is no good dehnition for aE from the holographic result 
because it may depend on the value of uj/m^ and as we will explain at the end of this 
subsection there will also be contributions to the DC conductivity at order cte from 
other quasinormal modes than the one that we focus here. 

There are several comments on this result which we list as follows. 


• The calculation above is valid for m? 3> w/ro, i.e. 1. If we want to 

go to the opposite limit <C tn/ro 1 in this holographic framework we can 
still perform this calculation except that the far region equations now are totally 
different. This can be seen from the second equation (4.2): now the coefficient 
in front of is more important. The result in this limit would recover our old 
result in [ 8 ], including the different form of cx^;. Thus one important conclusion is 
that the DC calculation in the main text based on radially conserved quantities 
is only consistent with the AC result in the limit ca/ro -C w?. 

• Note that here hi does not give us the value of because this is the <C 1 limit 
while Ts should be determined from the pole in the ca/ro ~ limit. In this limit, 
we can also solve the equations by assuming m? = Xiu/tq where Ai is an order 1 
number. In this limit, there will be order 0{u:) terms in the far region equations 
and we should solve these equations order by order in ca, i.e. 0(u) order solutions 
come from two parts: u order corrections to the 0{1) order solutions and one 
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linearly independent part. This is beyond the calcnlation in this paper and we 
will leave it for fntnre investigation. 

• The hydrodynamic formulas only capture the physics of the one quasinormal mode 
that we focus on, which has t^T —)• oo. However, the DC holographic result has 
contributions from all quasinormal modes. This means there will be an order 
1 /tsT difference in the holographic and hydrodynamic results for the DC conduc¬ 
tivity. 

• Finally note that these considerations go through with little change for the explicit 
breaking case. The result is given by (4.14) upon substituting = 2g^0Q. 


5 Conclusion and discussion 

We have considered two holographic models to encode axial charge dissipation in the 
probe limit. They are dual to four dimensional strongly coupled anomalous systems in 
presence of background magnetic held. In our hrst holographic model,the axial charge 
dissipation is realised by a charged scalar non-normalisable mode. At weak coupling 
this corresponds to introducing a fermionic mass term. In our second holographic 
model, the so called Stiickelberg massive U{1)a model, the U{1)a is broken by giving 
the axial current an anomalous dimension. As we have argued both models are closely 
related. Indeed if one choses the scalar mass = 0 and /i = ps = 0 the relevant 
equations coincide in both models. 

We found that in both these two models, positive magnetoconductivity is exactly 
quadratic in in the magnetic held strength. Moreover this remains true even in the 
case of small relaxation times when the axial charge can not be considered to be 
approximately conserved. 

This is consistent with the recent experiments [14, 15, 16, 17]. A recent weakly 
coupled theoretical proposal for ionic scattering [36] also found exact quadratic scaling 
of the magnetoconductivity with magnetic held. We note that our results, beyond 
being valid at strong coupling, are also diherent in that we consider high temperature 
and low chemical potentials. In particular our results suggest that at weak coupling 
but high temperature there should still be positive magnetoconductivity quadratic in 
B even when the fermi energy does not intersect any Landau level but lies in the gap 
region. Our models should also be of relevance for application to non-central heavy ion 
collisions quark gluon plasma where strong magnetic helds are present. 

We have introduced the axial charge relaxation by switching on constant sources. 
Anther interesting way to induce it should be via a random source that averages to 
zero, i.e. introducing disorder. We leave this possibility for future investigation. 
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A Appendix 

In this appendix we will list the details of the previons holographic resnlt on the con¬ 
ductivity for the case without axial charge dissipation [8] and the equations mentioned 
in the main text for reference. 


A.l The quantum critical conductivity cte in the holographic 
model without axial charge dissipation 


In this subsection, we recall the result of the quantum critical conductivity and its 
plot for reference [8]. From holography in the f/(l)y x 17(1)^ model with a background 
B field, we have 


(T = a E 


i {Waf 

Udp5/dfi5 


(A.l) 


where 


o'e 


dp5 

dp5 


vr^T 


, (3^ sec r-V^) , 

8 v2 / p |- 5 -a/i-/ 3^ j p |- 5 +a/i-^^ j 

^2^2 r [ 3-V^ j r [ 3+\/l^ 1 

A r 
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(A.2) 


(A.3) 


with [3 = The plot for the real part (cx^;) and the imaginary part of a can be 

found in Fig. 12. For large B, we have —)■ 0. 
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Figure 12: The real part and the imaginary part of a in the holographic model without 
axial charge dissipation [8]. 


A.2 Equations for the background in the explicit U{1)a break¬ 
ing model 


Substituting the ansatz (2.12) into (2.3-2.5) we find the the following background equa¬ 
tions of motion 


A': + w A, + 


J.2f 

Ml f'\ ! 2g^0^ 




V!' + -V! + ^ A 




r-" 
SBa 


f 


+ 


5 /' 

- + — 
r f 


+ 




q‘^Al m 


r'^P r^f r^f 


0 , 

0 , 

0 , 

0 , 

0 . 


(A.4) 

(A.5) 

(A.6) 

(A.7) 

(A.8) 


We focus on the solution which is regular on the black hole event horizon. Thus 
At{r = ro) = 0. We set V)(r = tq) = 0. Equations (A.6) and (A.7) can be simplified 
as {PVl + SBaAz)' = 0 and {PfV^ + SBaAt)' = 0. Thus we can integrate the two 
equations to get -|- SBaA^ = Cq and PfV^ + SBaAt = 0. Eq. (A.4) can be written 
as 


W' -p -A' _ — 


2 ,2 (85a) 

^ 


A = 0. 


(A.9) 


Note that I 4 can be totally fixed by adding normalisable boundary condition. 
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A.3 Equations for the fluctuations in the explicit U{1)a break¬ 
ing model 

A.3.1 Equations of motion for longitudinal fluctuations at zero momentum 


We have the following equations of motion for the longitunial fluctuations a*, a^, Uj, Uz, 0i, 02 
at zero momentum on the top of the background (2.12) 


^ ( 2gA0i + ^ 


+ (0 + y) ^ + 


8Ba , Aq^A^cj) 

V, 


r^f 
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U}a[ + + 2iqr‘^f[ - 020 ' + 00 ' 2 ) 
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J2 


Vz + 


v: + 


8Ba , 
8Ba 


+ (“ + - q^fAfj 


r^P 




2q^(l){Atat - aJA^) + 2iquj(t)2At 


+ (^ + y) p2 + (uj^ + q^Al - TiPPf - q'^fA^j 


tqu 


r4/2 7-4y2 


{4>ctt + 2At0i) 


0 , 

0 , 

0 , 

0 , 

0 , 

0 , 

0 . 


A.3.2 Equations of motion for transverse fluctuations at zero momentum 


The equations of motion for the transverse fluctuations ax, ay, Vx, Vy at zero momentum 
on the top of background (2.12) are the follows 




+ ( r + ^ ) + ( Z97 “ 2g^0^ ) j. ( Vpy + A'^ay 
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< + (^ f- 2 .V^) ^ (k'.. T Ay.) 
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0, (A.IO) 
0 (A.ll) 
0, (A.12) 
0 (A.13) 


After defining a± = a. ± iay, = n. ± iVy, we have 
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) =0, (A.14) 
1=0, (A.15) 


For zero density case Vt = At = Vz = = 0, hy repeating the calcnlation in the 

appendix of [ 8 ], we have a^x = o'yy = ttT and the Hall condnctivity axy = 0. When 
At = Vz = 0, i.e. with /is = 0 and /i 7 ^ 0, we have axx = o'yy = ttT and axy = 
which is the same as the case withont axial charge dissipation [ 8 ]. We do not have 
analytical solntions for other cases. 

A.3.3 Equations for DC conductivity calculation 

We have seven ODEs for the helds aronnd the backgronnd (2.12) in which we assnme 
the most general case with backgronnd /i and /is. The eqnations are the following 



8BaE 


2q‘^(j)‘^ar + 2g( - 020' + 002) = 0 


2q‘^(j)‘^ar + 2g( - 020' + 00'2) 




qttrcj) — q(t)a'^ — 2qar4>' = 0 
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A.3.4 Equations for transverse fluctuations at flnite u and k 


We consider transverse fluctuations Qx, ( 
on top of ( 2 . 12 ) and the equations are 
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A.4 Equations of fluctuations for massive U{1)a 


The linearised equations of motion for the fluctuations at general uj and k in the massive 
model read 
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(A.18) 

(A.19) 

(A. 20 ) 


with two constraint equations 
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